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Di-hadron correlations and jet-hadron correlations are frequently used to study interactions of 
the quark gluon plasma with the hard partons that form jets. The existing background subtraction 
methods for these studies depend on several assumptions and independent measurements of the 
Fourier coefficients of the combinatorial background. In this paper, we present a method for deter¬ 
mining the background using a fit to the reaction plane dependence of the background-dominated 
region of the near-side to extract the background. We also fit the of the background-dominated 
region of the near-side without the reaction plane dependence. To test the accuracy of these meth¬ 
ods, a simple model is used to simulate di-hadron and jet-hadron correlations with a combinatorial 
background similar to that observed in the data. The true signal is compared to the extracted sig¬ 
nal. The results are compared to results from two variants of the zero-yield-at-minimum (ZYAM) 
method. We test these methods for mid-peripheral and central collisions and for di-hadron and 
jet-hadron correlations. These methods are more precise than the ZYAM method with fewer as¬ 
sumptions about the shape and level of the combinatorial background, even in central collisions 
where the experimental resolution on the measurement of the reaction plane dependence is poor. 

These methods will allow more accurate studies of modifications of the away-side jet and will be 
particularly useful for studies of jet-hadron correlations, where the combinatorial background is 
poorly constrained from previous studies. 
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I. INTRODUCTION 

The quark gluon plasma (QGP), a strongly interacting 
liquid of quarks and gluons, is produced in high-energy 
nuclear collisions m- Hard probes such as jets are fre¬ 
quently used to study the QGP because they are pro¬ 
duced by hard scatterings early in the collision and prop¬ 
agate through the medium. Hard partons interact with 
the medium and lose energy, a process called jet quench¬ 
ing. 

The interactions of jets with the medium are com¬ 
monly studied using three experimental methods: mea¬ 
surements of single-particle spectra at highpx; di-hadron 
correlations where at least one particle is at high mo¬ 
mentum, and fully reconstructed jets. Observations of 
jet quenching at RHIG was one of the key signatures of 
the formation of the QGP. The initial observation re¬ 
lied on measurements of the nuclear modification factor 
RaAi which compares the hadron spectrum in A-\-A to 
that in p+p. At high momentum where hadron produc¬ 
tion is expected to be dominated by jets (px > 5 GeV/c 
at RHIG), the number of hadrons observed in A+A colli¬ 
sions is roughly 1/5 that expected fromp-|-p collisions [5l- 

m- 

Studies of jets in a QGP are complicated by the large 
background due to soft processes. Not only is there a 
large background but the strong collective flow observed 
in bulk particle production leads to correlations between 
particles in the bulk similar to the correlations due to 
jet production. Since both jets and collective flow con¬ 
tribute to correlations between particles, collective flow 


generates a significant background for any study of jets 
in heavy ion collisions. Gollective flow is dominant at low 
momenta (px 2 GeV/c), so this background has typi¬ 
cally been dealt with by focusing studies of jets on high 
momentum particles. However, gluon bremsstrahlung 
leads to gluons that are softer than the parent parton. 
As these hadronize, the final-state hadrons are softer on 
average than the final-state hadrons from the parent par- 
ton [To]. This means that many of these modifications 
are likely to be concentrated at low momentum and at 
large angles from the parent parton. Therefore, a reli¬ 
able and precise method for background subtraction is 
needed in order to quantify jet modihcation at low and 
intermediate momenta. 

We extract the background using a fit to the reac¬ 
tion plane dependence of background-dominated region 
on the near-side of di-hadron correlations and jet-hadron 
correlations, called the reaction plane fit (RPF) method. 
We also compare to a ht without the reaction plane fit, 
the near-side fit (NSF) method. We demonstrate that 
these methods produce more accurate and reliable re¬ 
sults than the ZYAM method. These methods both take 
advantage of differences between the signal and the back¬ 
ground in order to determine the background. The signal 
on the near-side is concentrated in a peak near the trig¬ 
ger hadron or jet. This peak is narrow in both azimuth 
and pseudorapidity. In contrast, the background forms 
a peak in azimuth but not pseudorapidity. For analy¬ 
ses in a narrow pseudorapidity range, the background 
is roughly independent of pseudorapidity. We combine a 
background roughly matching that observed in A+A col- 
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lisions with a known signal generated from PYTHIA m 
to show that these methods accurately and reliably re¬ 
construct the signal. 

All background subtraction methods for di-hadron cor¬ 
relations and jet-hadron correlations to date have as¬ 
sumed that the shape of the background is known. Our 
methods also make this assumption. Previous methods 
further assume that the magnitude of the coefficients of 
the Fourier decomposition of the background are known 
from other studies and can be measured independent of 
the correlation. Our methods do not make this assump¬ 
tion; effects that could modify these coefficients such as 
contributions from jets or differences in the hydrodynam- 
ical flow in events that contain jets, are taken into ac¬ 
count because the Fourier coefficients are fit. Moreover, 
this allows the accurate determination of the background 
even in cases where the Fourier coefhcients have not been 
measured to higher order, such as jet-hadron correlations. 
Like previous methods, our methods are also dependent 
on the assumption that contributions from other corre¬ 
lations such as Hanbury-Brown-Twiss correlations or de¬ 
cays of resonances that are not part of a jet are negligible. 
The ZYAM method makes the assumption that there is 
an angle in azimuth for which there are no correlations 
from jets. At sufficiently low momenta, this is certainly 
not true, since the near-side and away-side peaks overlap 
in azimuth. Our methods do not assume that the sig¬ 
nal is zero at a given angle, however, we assume instead 
that the signal is negligible on the near-side when the 
difference between the pseudorapidities of the associated 
particle and the trigger is large. 

We first summarize correlation studies, discussing pre¬ 
vious studies, correlations that contribute to the back¬ 
ground, the shape of the background, and the shape of 
the signal. This is a motivation for our model of the back¬ 
ground and the signal, discussed in the following section. 
We then show the results of the NSF method. While the 
NSF method is more accurate than the ZYAM method, 
the results are not stable when the fit range is reduced. 
The RPF method is tested for di-hadron correlations in 
mid-peripheral collisions with a fit over a wide and a nar¬ 
row range in azimuth, demonstrating that this method is 
more robust than the NSF method and produces more 
precise results than the ZYAM method. We then test 
this method for central collisions. Even though the reac¬ 
tion plane resolution is poor in central collisions, the lim¬ 
ited information available constrains the background and 
produces more precise results than the ZYAM method. 
Finally, the method is tested for jet-hadron correlations. 


II. CORRELATION STUDIES 

In a typical di-hadron correlation study [T^HTB] , a high- 
Pt trigger particle is identified and used to define the 
origin in azimuth and pseudorapidity. Typically, it is de¬ 
fined by its high momentum alone, restricted to a range 
of momenta. Here all trigger particles in a given mo¬ 


mentum region are accepted and then the correlation be¬ 
tween particles is studied in both azimuth and pseudo¬ 
rapidity. By selecting high-px particles the fraction of 
trigger particles coming from the production of jets is 
enhanced, however, it is not possible to determine con¬ 
clusively which trigger particles originate from hard pro¬ 
cesses and which arise from soft processes, particularly 
for lower-momentum (<10 GeV/c) trigger particles. 

Associated particles are also usually defined only as 
particles within a given momentum interval. For each 
associated particle in the event, its position relative to 
the trigger particle in azimuth (A0) and pseudorapidity 
(A? 7 ) is determined and the conditional yield is calcu¬ 
lated. In this paper, the conditional yield is normalized 
by the number of events. As with the trigger particle, it 
is not possible to determine conclusively which associated 
particles are from hard processes and which are from soft 
processes. For jet-hadron correlations, instead of a trig¬ 
ger hadron, a jet candidate is used to define the origin 
in azimuth and pseudorapidity HZIIIH]. Below, “trigger” 
refers to either a trigger particle or a trigger jet. 

Since both the associated particles and the trigger in¬ 
clude particles created from and modified by soft pro¬ 
cesses such as hydrodynamical flow, the combinatorial 
background is not azimuthally isotropic in heavy ion col¬ 
lisions. The way that this combinatorial background is 
typically treated is to assume that the contribution to 
both the trigger and associated particles can be factorized 
into a contribution from hard processes, the signal (J), 
and a contribution from soft processes, the background 
(B). This is referred to as the two-source model [Tn]. The 
term “raw signal” is used below for what would be mea¬ 
sured experimentally after corrections for detector and 
acceptance effects but before background subtraction. 
The raw signal contains signal-signal (J-J) correlations, 
signal-background [J-B) correlations, background-signal 
{B-J) correlations, and background-background {B-B) 
correlations. It is assumed that the processes that pro¬ 
duce the signal and the background are completely inde¬ 
pendent so that the J-B and B-J correlations are also 
background. 

A typical raw signal from di-hadron correlations for 
trigger momenta 8 < < 10 GeV/c within pseudora¬ 

pidities \r]\ < 0.5 and associated particles within |p| < 0.9 
with momenta 1.0 < < 2.0 GeV/c in p-\-p collisions 

at ^/s = 2.76 TeV in PYTHIA [TT] is shown in Figure]^ 
The raw rsignal is normalized by the number of equiva¬ 
lent Pb -j- Pb collisions because this is used as the known 
signal later and added to the background in Pb -\- Pb 
collisions. Because even PYTHIA has background from 
an underlying event, this raw signal includes J-J, J-H, 
H-J, and B-B correlations. Because there is no physi¬ 
cal correlation between the signal and the background in 
PYTHIA, the J-B, B-J, and B-B correlations are inde¬ 
pendent of azimuth and lead to the plateau in Figure 
Figure]^ shows that there is a peak near 0° which is nar¬ 
row in both A(j) and Arj. There is also a peak near 180°, 
which is narrow only in Acj), however, this peak is roughly 
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independent of pseudorapidity. The former is called the 
near-side and comes from associated particles from the 
same parton as the one that generated the trigger par¬ 
ticle. The latter is called the away-side and comes from 
associated particles from the parton that scattered off of 
the parton that generated the trigger particle. 

The parton that produces the near-side is generally 
thought to be biased towards partons that have not in¬ 
teracted strongly with the medium and it is therefore 
critical to study the away-side peak. Even in PYTHIA 
the away-side peak is roughly independent of pseudora¬ 
pidity within the typical acceptance used in correlation 
analyses, as shown in Figure While a hard parton 
scattering produces two back-to-back partons in the rest 
frame of the parton, the rest frame of the parton is, in 
general, not the same as the rest frame of the incoming 
nuclei. The difference in azimuth is negligible since most 
of the momenta of both the parton and the nuclei are in 
the direction of the beam pipe, however, the difference in 
pseudorapidity can be quite substantial. This causes the 
away-side to be broad in Arj without modified fragmen¬ 
tation or interaction with the medium. This is evident 
in Figure [2 

In a heavy ion collision, these peaks may be widened 
through partonic interactions with the medium |20j . 
for instance if the original partons have emitted 
bremsstrahlung gluons. The yield in the peaks may either 
be higher or lower, depending on how the parton inter¬ 
acted with the medium and on the specific momentum 
range. A parton that emitted a bremsstrahlung gluon 
would have less energy when it fragments, meaning that 
the peak at high associated particle momentum would 
be depleted because the odds of producing a high mo¬ 
mentum particle through fragmentation would be lower. 
However, at lower momentum, the peak would include 
particles from fragmentation of both the parent parton 
and the bremsstrahlung gluon, so the peak would be en¬ 
hanced at low associated momenta. 

For studies of the near-side, the raw signal at large 
Ar] can be used to determine the level of the background 
without extracting the precise [iiiEiiEg. However, 
this is only useful for studies of the near-side because the 
signal on the away-side is also roughly independent of Arj 
and will also be subtracted when this method is applied. 

The background due to soft processes can be written 
in general by a Fourier decomposition of the azimuthal 
anisotropy relative to the reaction plane: 


dN 

d{(j) - tpu) 


OO 

oc 1 -h ^ 2u„ cos[n{(l) - 'iPh)], 

n—1 


( 1 ) 


where N is the number of particles, (j) is the angle of a par¬ 
ticle’s momentum in azimuth in detector coordinates and 
ipH is the angle of the reaction plane in detector coordi¬ 
nates. In high energy heavy ion collisions the Fourier co¬ 
efficients Vn arise due to hydrodynamical flow [THl] . The 
initial overlap region is azimuthally anisotropic, leading 
to anisotropic pressure gradients, which give rise to the 



FIG. 1: Di-hadron correlations for trigger momenta 8 < py < 
10 GeV/c within pseudorapidities Ipl < 0.5 and associated 
particles within I p| < 0.9 with momenta 1.0 < <2.0GeV/c 

in p-fp collisions at y/s = 2.76 TeV in PYTHIA [TT|. The 
signal is normalized by the number of equivalent Ph + Pb 
collisions in our simulations and corrected for the acceptance 
using the mixed event correction described in the text. 


Vn- These initial azimuthal anisotropies are preserved 
through partonic hydrodynamical flow and lead to az¬ 
imuthal anisotropies in the final-state hadrons. The mag¬ 
nitude of the Fourier coefficients decreases with in¬ 
creasing order. The sign of the first-order coefficient vi 
is dependent on the incoming direction of the nuclei and 
changes sign when going from positive to negative pseu¬ 
dorapidities. Since correlation analyses typically average 
over both positive and negative pseudorapidities, the av¬ 
erage vi is zero. 

The even Vn are generally understood to arise mainly 
from anisotropies in the average overlap region of the in¬ 
coming nuclei, considering the nucleons to be smoothly 
distributed in the nucleus with the density depending 
only on the radius. The with even n are correlated 
with the reaction plane. The odd are generally un¬ 
derstood to arise from the fluctuations in the positions 
of the nucleons within the nucleus. High-energy heavy 
ion collisions happen on a time scale short enough to be 
sensitive to the position of individual nucleons within the 
nucleus. Since these fluctuations are not causally related 
to the reaction plane, the odd are assumed to be un¬ 
correlated with the reaction plane. Recent measurements 
by ATLAS confirm that the correlation between n = 2 
and n = 3 reaction planes is very weak |23j . 

For B-B correlations entirely due to hydrodynamical 
flow the conditional yield will be given by [24] : 

7 

cos(nA</))], (2) 

where H is a constant that depends on the multiplic¬ 
ity of trigger and associated particles in an event and 
on the normalization convention. At/) is the difference 
in azimuthal angle between the associated particle and 
the trigger, Vn is the v„ for the trigger, and Vn is the 
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FIG. 2: Di-hadron correlation signal for trigger momenta 
8 < py < 10 GeV/c within pseudorapidities \ri\ < 0.5 and as¬ 
sociated particles within |» 7 | < 0.9 with momenta 1.0 < < 

2.0 GeV/c in 30-40% Pb + Pb collisions at = 2.76 TeV. 

The signal is normalized by the nnmber of Pb + Pb collisions. 
The signal is from Figure and the generation of the back¬ 
ground is discussed in the method section. 


Vn for the associated particle. B-B correlations due to 
processes other than hydrodynamical flow are generally 
assumed to be negligible. In a typical analysis, the pseu¬ 
dorapidity range for both trigger and associated particles 
is restricted to a region where the do not change dra¬ 
matically within the acceptance for the analysis and in 
this case the pseudorapidity dependence of ^ is neg¬ 
ligible. We consider only such analyses here, although 
analyses over a wide enough range in pseudorapidity for 
the Vn to change significantly are possible m- The shape 
of a typical di-hadron correlation for trigger momenta 8 

< Py < 10 GeV/c within pseudorapidities \r]\ < 0.5 and 
associated particles within |p| < 0.9 with momenta 1.0 

< < 2.0 GeV/c in 30-40% Pb + Pb collisions at ^snn 
= 2.76 TeV including background terms up to n = 
10 is shown in Figure The signal is normalized by the 
number of Pb + Pb collisions. The signal is from Figure[2 
and the generation of the background is discussed in the 
Sec. III. 

In collision systems with no hydrodynamical flow 
where the background is due to the underlying event or 
soft processes uncorrelated with the reaction plane, the 
J-B correlations will be independent of At/) and there¬ 
fore only add an overall constant background. This as¬ 
sumption has been used for studies of di-hadron corre¬ 
lations in p+p and d + Au collisions [m where 

this constant term is assumed to arise from the underly¬ 
ing event. However, jets are correlated with the reaction 
plane because jets are quenched more out-of-plane, where 
the mean path length of medium a parton must traverse 
is longer, than in-plane p5H5T| . Therefore, when there 
is a background due to hydrodynamical flow and a jet 
signal suppressed by jet quenching, the signal and the 
background are both correlated with the reaction plane. 
The J-B and B-J correlations will not be independent 
of A(/) in this case. Since it is always possible to write 


any function as a Fourier decomposition. Equation also 
describes J-B and B-J correlations. This leads to an 
overall background due to J-B, B-J, and B-B correla¬ 
tions given by 

7 

—— = b[1 + Y^ cos(nA(/))], (3) 

^ n—1 

where Vn {Vn) is the pair weighted average of the Vn (Vn) 
due to jet quenching and the z;* (n^) due to hydrody¬ 
namical flow. 

Generally the used in background subtraction are 
measured separately from the correlation measurements. 
The appropriate method for measuring the is not ob¬ 
vious. Different methods for measuring flow produce sys¬ 
tematically different results in the same event class. Mea¬ 
surements of Vn using the event plane lead to systemati¬ 
cally higher results than measurements using correlations 
between multiple particles, such as a four-particle cumu- 
lant method [32]. The latter is less sensitive to non-flow, 
making it less sensitive to contamination from jets. In 
principle this would make it a better measurement for the 
Vn in Equation and Equation however, these meth¬ 
ods are also less sensitive to event-by-event fluctuations 
in flow and to local hot or cold spots in the medium. Ad¬ 
ditionally, it is possible that events containing jets could 
have slightly different average Vn due to hydrodynamical 
flow than measurements of Vn in minimum bias collisions. 
For these reasons it is desirable to have a method where 
the Vn and Vn used for the background are determined 
from the same analysis as the measurement of the J-J 
correlations. 

Furthermore the determination of the appropriate B in 
Equation is difficult and prone to assumptions about 
the signal which may not be true. The most common 
method used is to assume a zero yield at minimum 
(ZYAM) [23|33H35], or some variation of ZYAM, for in¬ 
stance to assume zero yield at A(/) = 1. This assumes 
that there is a region in A^ where the signal goes to 
zero. The problem with this assumption is that there 
may be no region in azimuth where the J-J correlations 
go to zero. Even in PYTHIA, at low momentum < 
1 GeV/c) there is no flat region in Ai/), indicating that 
there is no reliable At/) region where the signal can be 
assumed to be zero. In heavy ion collisions, where both 
the near-side m and the away-side |3S] peak may be 
modified by interactions with the medium, it is even less 
reliable. The ABS method [37] uses mixed events to de¬ 
termine the background level. This is an improvement 
on ZYAM, however, it makes the assumption that the 
number of J-B and B-J pairs are negligible compared to 
the number of B-B pairs. This assumption is valid for 
central collisions where the multiplicity of background 
particles is large, but not valid for peripheral collisions 
or collisions in small systems. 

The data at small A^ and large Ap are observed to 
be dominated by the background II1I1I1I12], whatever 
its source, and can be fit to Equation to determine 
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the v^vl^. The background determined in this manner 
still assumes that the form of the background in Equa¬ 
tion is correct and it is sensitive to the validity of the 
assumption that there is no residual signal in the large 
A 77 and small A(/) region used for these fits. However, 
it improves on ZYAM because there is no assumption 
that the signal goes to zero at one point and improves 
on both ZYAM and the ABS method because there is no 
assumption that the in the background in correla¬ 
tion studies are equal to the measured from studies 
of hydrodynamical flow. 


III. METHOD 

We focus on di-hadron correlations and jet-hadron cor¬ 
relations with associated particle momenta 1 < < 

2 GeV/c within |ry| < 0.9 for both 0-10% and 30-40% 
central Pb + Pb events at i/snn = 2.76 TeV. We select 
trigger hadrons with 8 < < 10 GeV/c with \ri\ < 

0.5 and trigger partons with 20 < p^ < 40 GeV/c with 
\ri\ < 0.5. While the ALIGE detector can select trigger 
hadrons over a wider 77 acceptance, using the same 77 se¬ 
lection for trigger hadrons and partons simplified the sim¬ 
ulations of the background. The signal is generated from 
PYTHIA [TT] events using the Perugia 2011 tune [38] . 

Di-hadron correlations are calculated using charged 
hadrons for both the trigger and associated particles. 
Jet-hadron correlations are calculated using gluons and 
quarks as a proxy for fully reconstructed jets. We do not 
attempt a realistic simulation of jets. While real data 
could lead to fake jets, they would comprise particles 
correlated by flow and therefore would not have an as¬ 
sociated near-side peak. Instead, fake jets would change 
the hjj. Since we fit the h*, we would extract the correct 
for background subtraction using our fit. 

The background is generated assuming that each trig¬ 
ger and each associated particle is correlated with the 
reaction plane with the Vn given in Table |l] up to = 
10. The available data guide the choice of Vn [M1I39H4T] . 
The exact choice of Vn does not impact whether or not 
the method is feasible. Larger Vn make the background 
more difficult to extract, particularly for the higher-order 
Vn where the available data do not constrain the Vn sig¬ 
nificantly. We therefore use upper bounds in order to 
test the method in a worst case scenario. Eor associ¬ 
ated particles and trigger particles for di-hadron correla¬ 
tions, 772 , 773 , 774 , and 775 are chosen to approximate the 
values observed in the data gni EH- The data avail¬ 
able do not tightly constrain the higher-order Vn so we 
use 77 ji +2 = 77 „/ 2 . This is an approximate upper bound. 
Only 772 is available for reconstructed jets m- We esti¬ 
mate that the 77 „ for jets is approximately the same as 
the Vn for high-px hadrons. 

To get the signal to background correct, 1660 PYTHIA 
events are simulated for each 0-10% central Pb+Pb event 
and 251 PYTHIA events for each 30-40% central Pb+Pb 
event, the number of binary nucleon-nucleon collisions in 


TABLE I: 77,1 values used for background calculations. For 
77n>5 we use 7;„+2 = v„/2. 



V2 

773 

774 

775 

0-10% 

Vn 

vi 

0.041 

0.030 

0.030 

0.030 

0.0023 

0.0150 

0.0011 

0.0100 

30-40% 

Vn 

0.134 

0.047 

0.0173 

0.0092 

vi 

0.100 

0.030 

0.0150 

0.0100 


each Pb + Pb collisions calculated by CMS [H] . Since 
PYTHIA events include an underlying event, we subtract 
this background using the ZYAM method and use this 
as our known signal. We emulate approximately 8 x 10® 
0-10% and 16 x 10® 30-40% central Pb -I- Pb collisions. 

To simulate the background pairs, the reaction plane 
angle 'ip is chosen to be zero in detector coordinates and 
random trigger jets or hadrons and associated particles 
are thrown with a distribution described by Equation 
with the Vn given in Table |Tj We estimate the number of 
associated particles with 1 < < 2 GeV/c with I 77 I < 

0.9 using charged hadron gSl El spectra measured by 
ALICE. We throw a random 77 for both the associated 
particle and the trigger. To emulate the approximate 
effect of 77 dependent 77„, we apply a 1% linear decrease 
in the 77 „ from 77 = 0 to 77 = 0.9, consistent with the slight 
77 dependence observed at midrapidity. The observed Vn 
will be largest when the Trth-order event plane is used, 
however, in a typical analysis to measure the correlations 
due to jets, the second-order event plane is used. In our 
model we assume that all event planes for even n are 
identical to the second-order event plane. There is no 
correlation between the odd and even ?r-event planes. We 
therefore choose a random orientation for the odd 77-event 
plane for each simulated event. 

The RPF method uses the reaction plane dependence 
to determine the background. The finite resolution for 
reconstructing the reaction plane changes the shape of 
the background for reaction plane dependent correlation 
studies Plj. To simulate a realistic measurement, the 
true reaction plane angle ip is smeared with a Gaussian 
with a width of 20° for 30-40% central collisions and 40° 
for 0-10% central collisions. This reaction plane resolu¬ 
tion is quantified in terms of 

Vn =< COs[n{lptrue “ i’reco)] > ■ (4) 

For perfect reaction plane reconstruction r„ = 1 and for 
no reaction plane resolution r„ = 0. In our model we get 
r 2 = 0.79, r 4 = 0.38, and vq = 0.11 for 30-40% collisions 
and r 2 = 0.58, r 4 « 0, and re « 0 for 0-10% collisions. For 
odd n, Vn = 0 because the odd and even n reaction planes 
are not correlated. Data indicate that the even reaction 
planes are not 100% correlated [23]. This changes the 
effective 77 „ when an analysis is done for a trigger fixed 
relative to the reaction plane, however, this is taken into 
account by using the r„ measured relative to the reaction 
plane used in the analysis. 
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FIG. 3: Acceptance correction for a trigger with a flat distri¬ 
bution within \ri\ < 0.5 and an associated particle with a flat 
distribution within \rj\ < 0.9. 


The sharp cut off in the pseudorapidity r] of the par¬ 
ticles accepted leads to a trivial structure unrelated to 
physics. Pairs with Ary = ry* — ry“ « 0 have nearly 100% 
acceptance, however, pairs with |Ary| « rjmax'^Vmax have 
nearly 0% acceptance. In measurements of correlations, 
this is corrected with a mixed event correction which in¬ 
cludes detector effects. While our model has no detector 
effects, this trivial acceptance effect significantly mod¬ 
ifies the simulated signal. We therefore also apply an 
acceptance correction. If ry^^^, ^ Vmax^ there will be a 
plateau between - 17 y^„,^-ry^„,, I and \vLax-V^ax\- This 
is shown in Figure for a trigger with a flat distribution 
within Iryl < 0.5 and an associated particle with a flat 
distribution within | 7 y| < 0.9. 


IV. RESULTS 

We first test the NSF method. The raw signal in Fig¬ 
ure is projected from 1.0 < Ary < 1.4. This is then 
normalized by the Ary width of this projection in order to 
retain roughly the same scale as in Figure]^ independent 
of the range of the projection. The extracted background 
is compared to the true background and the background 
extracted using a variant of the ZYAM method. The sig¬ 
nal is then extracted using the fit background and com¬ 
pared to the true signal and two variants of the ZYAM 
method. 

The same procedure is followed for the reaction plane 
dependence, testing the method with different fit ranges, 
for 30-40% central collisions, and for di-hadron and jet- 
hadron correlations. The same Ayy range and normaliza¬ 
tions are used for the projections. In order to make the 
discussion easier to follow, we use di-hadron correlations 
in 30-40% central collisions with a fit range of |A0| < 
tt/2 as our primary reference and only vary one condi¬ 
tion for the fit at a time. For each sample, the true and 
extracted backgrounds are compared and then the true 
and extracted signals are compared. The same symbols 


are used throughout the discussion for clarity. 

The four methods used for the background subtraction 
are: 

• The ZYAl method, a variation of ZYAM where the 
background is fixed at A(j) = 1 instead of at the 
minimum; 

• The modified ZYAl method, a variation of ZYAl 
where the background is fixed using only data in 
the background-dominated region, 1.0 < |Ary | < 
1.4; 

• The NSF method, which fits the near-side in 
the background-dominated region to determine the 
background; 

• The RPF method, which fits the reaction-plane- 
dependent near-side in the background-dominated 
region to determine the background. 

The ZYAl method is less sensitive to statistical fluctu¬ 
ations than ZYAM. The modified ZYAl is less sensitive 
to the signal than ZYAl and the background determined 
from this method can be directly compared to the NSF 
and RPF methods because they all use the same data. 
The ZYAl method requires the as input. In an anal¬ 
ysis of data, the are typically taken from other stud¬ 
ies and the have an uncertainty. In our analysis we 
assume that this uncertainty is 5%, comparable with un¬ 
certainties on Vn measured with a single method, assume 
that the uncertainties for the trigger and associated par¬ 
ticles are correlated, and use the true value as the nom¬ 
inal value. Since the methods for measuring vary in 
their sensitivity to fluctuations and non-flow contribu¬ 
tions, this likely underestimates the uncertainty on the 
proper to use for the background in a di-hadron cor¬ 
relation measurement. In addition, the nominal value of 
Vn used for background subtraction is not likely to be 
centered at the exact true Vn- This therefore likely un¬ 
derestimates the uncertainties in the ZYAl method. 

We compare the true yield to the yield extracted using 
various methods. The yield is given by 

dN /•'’ d^N 
^ ~ NedAr] ~ X AedA^dA^^^*^ 

where a = -1.05 and b = 1.05 for the near-side and a 
= 2.09 and b = 4.19 for the away-side. The yields are 
not comparable between 0-10% and 30-40% central data 
because of the normalization by the number of events. 

A. The Near-Side Fit Method 

The raw signal in the region 1.0 < \At] \ < 1.4 is fit to 
Equation]^ to order n = 3 from |A(/)| < 7 r/ 2 . Figure]^ 
shows the true background, the signal plus background, 
and the background from this fit for di-hadron correla¬ 
tions in 30-40% Pb + Pb collisions at Asnn = 2.76 TeV. 
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FIG. 4: Top: Signal + background for di-hadron correlations 
in 30-40% Pb+Pb collisions at .^snn = 2.76 TeV in the region 
1.0 < IA 77 I < 1.4. This is compared to the true background, 
the background from the modified ZYAl method, and the 
background from the NSF method. (See text for details.) 
The ht for the NSF method is to Equation to order n = 
4 from |A<()| < n/2 and has x^/NDF = 63.6/45. Bottom: 
Ratios of the background from the NSF and ZYAl methods 
to the true background. 


FIG. 5: Top: The true signal for di-hadron correlations in 
30-40% Pb + Pb collisions at = 2.76 TeV. This is com¬ 

pared to the signal extracted using the background from the 
ZYAl method, the modified ZYAl method, and the back¬ 
ground extracted from the NSF method for | At/I < 7 r /2 using 
the fit shown in Figure (See text for details.) Bottom: 
Differences between the true signal and the signal extracted 
using the background from the ZYAl method, modihed ZYAl 
method, and the background from the RPF method. 


Figure]^ also shows the background extracted in this re¬ 
gion using the modified ZYAl method. 

The signal extracted using the NSF background, the 
modified ZYAl background, and the standard ZYAl 
background are compared in Figure to the true sig¬ 
nal. Only statistical uncertainties on the background are 
shown. Figure]^ and Figureshow that the NSF, ZYAl, 
and modified ZYAl methods describe the background 
well in this model and that they have comparable un¬ 
certainties. The nominal value of the signal extracted 
using the fit is not centered on the true value like the 
nominal values for the ZYAl and modified ZYAl meth¬ 
ods, however, as discussed above, our implementation of 
these methods may slightly underestimate the true uncer¬ 
tainties on the Vn and therefore underestimate the shape 
distortions from the ZYAl method. Additionally, in a 
heavy ion collision, the away-side could be modified sig¬ 
nificantly and become much broader. In this scenario, 
the NSF method would be less sensitive to residual sig¬ 
nal in the away-side than the ZYAl method. 

Figure and Figure compare the signal extracted 
using the ZYAl method, modified ZYAl method, and 
NSF method for \A(j)\ < 7 r /2 for di-hadron correlations 
in 0-10% central Pb + Pb collisions and jet-hadron cor¬ 
relations in 30-40% central Pb + Pb collisions. The NSF 
method describes the background better for 0 - 10 % cen¬ 
tral collisions than for 30-40% central collisions. In cen¬ 
tral collisions the u* are smaller and the background is 
larger so it is possible to determine the background with 


higher precision. The yields are extracted for all these 
cases using Equation and are summarized in the Ta¬ 
ble [n] The NSF method provides a more precise mea¬ 
surement of the yield than ZYAl in all cases. 

However, we noticed that the fit was sensitive to 
the fit range. Figure shows the true background, 
the signal-{-background, the background using the ZYAl 
method, and the background extracted using the NSF 
method in the range |A(/)| < 1.25. This fit clearly fails 
to describe the background, even though the fit con¬ 
verged and the fit quality is comparable to the fit in the 
range | A(()| < tt/2. This is because the fit needs to be able 
to distinguish between V 2 and V 3 . At Acj) = tt/S, the V 3 
term starts increasing while the V 2 term is still decreas¬ 
ing. Either the fit needs to cover enough range in At/) to 
distinguish between these two terms or the data need to 
have enough statistics that it is possible to discern the 
relative weights of the V 2 and U 3 terms from the width 
of the near-side peak. The fit in Figure shows that 
realistic statistics do not provide data with the precision 
required for the latter. If a fit to the near-side at large At] 
were used for analyzing data, it could lead to subtracting 
the wrong background and potentially extracting a sig¬ 
nal with a distorted signal on the away-side. We there¬ 
fore explore using the reaction plane dependence of the 
raw signal, which uses more information and therefore is 
likely to have lower uncertainties and be more stable. 








FIG. 6: Top: The true signal for di-hadron correlations in 0- 
10% Pb + Pb collisions at y%NN = 2.76 TeV. This is compared 
to the signal extracted using the background from the ZYAl 
method, the modified ZYAl method, and the background ex¬ 
tracted from the NSF method for |A(;li| < tx/2. Bottom: Dif¬ 
ferences between the true signal and the signal extracted us¬ 
ing the background from the ZYAl method, modified ZYAl 
method, and the background from the RPF method. 


TABLE II: Yields Y as defined in Equation [^scaled by 10“^ 
from FigureFigureand Figure]^ The first uncertainty 
on the true yield is the statistical uncertainty and the second 
is the uncertainty due to the uncertainty on the subtraction 
of the background from the underlying event. The first un¬ 
certainty for the ZYAl method is the statistical uncertainty, 
the second uncertainty is the uncertainty on the background 
level, and the third uncertainty is the uncertainty due to uf 
and V 2 - ZYAl uncertainties are propagated assuming 100% 
correlation between V 2 and V 2 and no correlation between the 
uncertainty on the level of the background and the uncertainty 
on the Vn- Errors due to higher order Vn are not considered 
but are approximately 10% of the uncertainties due to Sf a-nd 


Sample 


Yield (y 
near-side 

X 10“^) 

away-side 

30-40% 

h-h 

True 

Mod. ZYAl 
Std. ZYAl 
NSF 

17.1 ± 0.1 ± 0.2 
18.9 ± 4.2 ± 1.2 
15.7 ± 1.6 ± 1.2 
17.14 ± 1.1 

19.9 ± 0.1 ± 0.2 

21.9 ± 4.2 ± 1.2 
18.7 ± 1.6 ± 1.2 

20.14 ± 1.11 

0-10% 

h-h 

True 

Mod. ZYAl 
Std. ZYAl 
NSF 

114.4 ± 0.2 ± 0.8 
75.5 ± 18.3 ± 0.9 
86.7 ± 7.0 ± 0.9 
111.63 ± 3.01 

132.8 ± 0.2 ± 0.8 
95.7 ± 18.3 ± 0.9 

106.9 ± 7.0 ± 0.9 
131.82 ± 3.01 

30-40% 

jet-h 

True 

Mod. ZYAl 
Std. ZYAl 
NSF 

13.19 ± 0.04 ± 0.17 

16.2 ± 4.2 ± 1.2 

13.2 ± 1.6 ± 1.2 
13.13 ± 0.77 

4.96 ± 0.04 ± 0.17 

8.2 ± 4.2 ± 1.2 

5.2 ± 1.6 ± 1.2 
5.13 ± 0.78 


C0.005 


True signal 
NSF signal 

Standard ZYAl signal 
Modified ZYAl signal 



FIG. 7: Top: The true signal for jet-hadron correlations 
in 30-40% Pb -F Pb collisions at = 2.76 TeV. This 

is compared to the signal extracted using the background 
from the ZYAl method, the modified ZYAl method, and 
the background extracted from the NSF method for |A<j!>| < 
tt/2. Bottom: Differences between the true signal and the sig¬ 
nal extracted using the background from the ZYAl method, 
modified ZYAl method, and the background from the RPF 
method. 
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FIG. 8: Top: Signal+background for di-hadron correlations 
in 30-40% Pb+Pb collisions at ^snn = 2.76 TeV in the region 
1.0 < IA77 I < 1.4. This is compared to the true background, 
the background from the modified ZYAl method, and the 
background from the NSF method. (See text for details.) 
The fit for the NSF method is to Equation to order n = 
4 from |A<(>| < 1.25 and has x^/NDF = 50.8/35. Bottom: 
Ratios of the background from the NSF and ZYAl methods 
to the true background. 
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B. The Reaction Plane Fit Method 

When a trigger hadron or parton is restricted relative 
to the reaction plane, the level of the background and 
the effective v\ are affected. The derivation of the ap¬ 


propriate reaction plane dependent forms are discussed 
in Ref. [21]. When the trigger is restricted to a range of 
angles relative to the reconstructed reaction plane, the 
effective even are given by 


r,R,t — 


Vn + cos(n(^s)' 


i(nc) 


Rr, 




fc=2.4,6. 


^ {Vk+n + ) COs(fc(^s) ' 


n{kc) 

kc 


R„ 


1 + Efc=2.4,6... cos{k(j)s)^^^^R, 


( 6 ) 


and the effective background level is given by 

/3fl = l+ ^ 2vk cos(fc^s) Rn (7) 

fc=2,4,6... 

where (ps is the center of range and 2 c is the width of the 
range [21]. The background is then given by: 

7 /\r 

= pR{l + cos{nA(j))). (8) 

Since the reaction planes for odd n are uncorrelated with 
the n = 2 reaction plane, all odd n terms have = 
vl^ when the n = 2 reaction plane is used for an anal¬ 
ysis. Here we consider simultaneous measurements of 

W dA<i>dArj trigger restricted to four different re¬ 

gions: 

• All; 

• In-plane: (l)s = 0: c= 7 r/ 6 ; 

• Mid-plane: (j)s = 7 r /4 and (j)s = 37 r/ 4 , c= 7 r/ 12 ; 

• Out-of-plane: (j>s = 7 r/ 2 , c= 7 r/ 6 . 

Note that the mid-plane range is actually split into four 
symmetric regions. These regions are shown schemati¬ 
cally in Figure The information on the reaction plane 
dependence of the raw correlations can reduce the un¬ 
certainty on the background. This can be understood 
by considering what additional information the reaction 
plane dependence provides. The level of the in-plane cor¬ 
relation is increased by and the term is increased, 
as shown in Equation and Equation This allows V 2 
to be determined with high precision from these corre¬ 
lations. In contrast, the background level of the mid¬ 
plane correlation is insensitive to V 2 and the modulation 
of the correlation by V 2 ’ is approximately equal to 
Since the in-plane and out-of-plane correlations strongly 
constrain the mid-plane correlation can be used to 
constrain the higher-order Vn- The normalization of the 
correlations per event allows the constant j3ji to be the 
same for all reaction plane orientations. 



FIG. 9: Schematic diagram showing the reaction plane angles 
used in the analysis. 


1. Di-hadron correlations in 30-40% central Pb + Pb 
collisions with the RPF method over |A(;A| <t:/2 


Figure shows the signal-l-background in the region 
1.0 < |A ?7 I < 1.4 for di-hadron correlations in 30-40% 
Pb-\-Pb collisions at y%NN = 2.76 TeV for in-plane, mid¬ 
plane, out-of-plane, and all reaction plane angles. We fit 
the signal on the near-side to all reaction plane orienta¬ 
tions simultaneously, restricting the fits for each reaction 
plane orientation to |A0| < 7 r/ 2 . The are fixed at the 
values in our model. In an experimental analysis, the r„ 
can be measured. We varied the order of Vn used in the 
fit until we used the fewest parameters necessary to get 
a reliable fit. We found that the fit worked best to n = 
4, corresponding to a total of six parameters, B, V 2 , V 2 , 
v^, U 4 , and 04 . The extracted parameters were all within 
error of the parameters used in the simulation. The RPF 
method is compared to the modified ZYAl method. Note 
that for the ZYAl method a different v!^ must be used for 
each reaction plane orientation. Again we use the nomi¬ 
nal values of thrown using our model and assume 5% 
uncertainties on the u*. Only statistical uncertainties are 
shown for each background method. Figure [TT] compares 
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the signal extracted using the various background meth¬ 
ods to the true background. The yields are given in Ta¬ 
ble HU The RPF method leads to a much more accurate 
determination of the signal shape and the yield than the 
ZYAl or modified ZYAl methods for correlations with a 
trigger in-plane and out-of-plane. This is because even 
with the optimistic 5% uncertainty on the the amplifi¬ 
cation of the V 2 for these correlations makes the accurate 
determination of the background difficult. 


2. Di-hadron correlations in 30-40% central Pb + Ph 
collisions with the RPF method over |A<^| < 1 


One of the benefits of the simultaneous fit of the back¬ 
ground to the different reaction plane orientations on the 
near-side is that it may allow a fit over a narrower A(^ 
range, which would reduce the impact of any signal from 
the away-side from a modified away-side. We therefore 
fit the signal over |A0| < 1, which should be less sensi¬ 
tive to residual signal than a fit to |A(/)| < 7r/2. This is a 
narrower range than that shown in Figure because by 
varying the fit range we found that the simultaneous fit 
converged even when a narrower range in A(/) was used. 
The background extracted with this fit is compared in 
Figure to the true background and the background 
extracted using the modified ZYAl method. The sig¬ 
nal extracted using this fit is shown in Figure and 


the yields are given in Table III While the signal using 
the fit over |A(/)| < 7r/2 shown in Figure 11 has slightly 
smaller uncertainties, the signal using the fit over |A0| < 
1 shown in Figure is comparable. Since the width of 
the away-side peak varies with both and and even 
the near-side width is observed to be modified in heavy 
ion collisions HU, a narrower fit range is better. With a 
range covering |A(/)| < 1, even an away-side peak with a 
width of A(j) = 1 would have less than 5% of its amplitude 
in the fit region. 


3. Di-hadron correlations in 0-10% central Ph-\- Pb 
collisions with the RPF method over |A(^| <71-/2 

Central collisions are often considered the most inter¬ 
esting because the medium reaches higher energy densi¬ 
ties and hottest temperatures. Naively one might assume 
that this method cannot be applied to central collisions 
because the reaction plane is not known to better than 
approximately 40° and the reaction plane bins described 
above are 30° wide. Equationj^and Equationj^show that 
the /3^ and differ less between in-plane, mid-plane, 
and out-of-plane when the are small than when the 
r„ are large. This is because a large fraction of the trig¬ 
ger particles reconstructed in-plane will actually come 
from mid-plane and some will even come from out-of- 
plane. Still, there is a difference between the correlations 
for trigger particles reconstructed with different reaction 
plane angles and this provides some information. 


The background determined using the ZYAl method, 
modified ZYAl method, and RPF method over the range 
|A^| < •7r/2 are compared to the signal-bbackground and 
the true background in our model for di-hadron correla¬ 
tions in 0-10% Pb-\-Pb collisions at ^snn = 2.76 TeV for 
in-plane, mid-plane, out-of-plane, and all reaction plane 
angles in Figure The fit used in the RPF method is 
to order n = 3 because the fit to order n = 4 was not 
stable. The reaction plane dependence of the correlations 
provide information on ■U 2 , however, because r 4 « 0, there 
is no additional information on V 4 . Figure [TS] compares 
the signal extracted with the ZYAl, modified ZYAl, and 
RPF methods and the yields are given in Table |III| This 
shows that the RPF method determines the signal with 
much higher accuracy and precision than the ZYAl or 
modified ZYAl methods. 


4- Jet-hadron correlations in 30-40% central Pb-\- Pb 
collisions with the RPF method over |A(^| < 7t/2 


We also explore using the reaction plane dependence 
for determination of the background in jet-hadron corre¬ 
lations. Figure [T6| shows the signal-1-background in our 
model in the region 1.0 < jAry | < 1.4 for jet-hadron 
correlations in 30-40% Pb -\- Ph collisions at ^snn = 
2.76 TeV for in-plane, mid-plane, out-of-plane, and all 
reaction plane angles and compares the background from 
the RPF method and the background from the modified 
ZYAl method. Figure compares the true signal to 
the signal extracted from using the the RPF method, 
ZYAl, and modified ZYAl methods to determine the 


background and the yields are given in Table III The 
RPF method works best. This method would be par¬ 
ticularly useful for jet-hadron correlations because only 
V 2 has been measured for jets and therefore background 
subtraction requires either large estimates of the uncer¬ 
tainty due to the m or limiting the analysis to the 
near-side where the background can be determines from 
the Ar] dependence of the signal [IH] . 


V. CONCLUSIONS 

We have presented two methods for determining the 
combinatorial background in di-hadron correlations and 
jet-hadron correlations by fitting the raw correlation at 
large Ary where the correlation is background-dominated. 
We have demonstrated that these methods accurately 
and reliably subtract the background using a model 
where the background is entirely due to flow and the sig¬ 
nal is generated using PYTHIA. These methods produces 
better results than the ZYAl method. The RPF method 
is more accurate than ZYAl even when the A(j) range of 
the fit is restricted and even in central collisions where the 
reaction plane resolution is poor. These methods will be 
particularly useful for jet-hadron correlations since only 
V 2 has been measured for jets, limiting the application of 
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FIG. 10: Top: Signal-|-background for di-hadron correlations in 30-40% Pb + Pb collisions at = 2.76 TeV in the region 

1.0 < \Ari I < 1.4 for in-plane, mid-plane, and out-of-plane triggers and for all triggers combined. This is compared to the 
true background, the background from the modified ZYAl method, and the background from the RPF method. (See text for 
details.) The data for all angles relative to the reaction plane have been scaled by 1/3 in the top panel in order to fit on the 
same scale. The fit for the RPF method is to Equationj^to order n = 4 from | A(/>| < tt/2 and has x^/NDF = 176/138. Bottom: 
Ratios of the background from the RPF and ZYAl methods to the true background. 
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FIG. 11: Top: The true signal for di-hadron correlations in 30-40% Pb + Pb collisions at y%NN = 2.76 TeV for in-plane, mid¬ 
plane, and out-of-plane triggers and for all triggers combined. This is compared to the signal extracted using the background 
from the ZYAl method, the modified ZYAl method, and the background extracted from the RPF method for \A(j)\ < 'k/2 

(See text for details.) The data for all angles relative to the reaction plane have been scaled 
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using the fit shown in Figure 

by 1/3 in the top panel in order to fit on the same scale. Bottom: Differences between the true signal and the signal extracted 
using the background from the ZYAl method, modified ZYAl method, and the background from the RPF method. 


methods such as ZYAM, ZYAl, or the ABS method [37] 
which all require the as input. The reaction plane 
fit method makes the same assumptions about the shape 
of the background made in other methods, namely that 
it has the functional form given by Equation How¬ 
ever, it does not make assumptions about the level of the 


background B or the The primary assumptions of 
this method are that the background has the functional 
form given by Equation]^ that the background’s Arj de¬ 
pendence is negligible, and that the contribution of the 
signal to the correlation at large Arj and small A(j) is 
negligible. This latter assumption is not valid at lower 
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FIG. 12: Top: Signal-|-background for di-hadron correlations in 30-40% Pb + Pb collisions at = 2.76 TeV in the region 

1.0 < \Ari I < 1.4 for in-plane, mid-plane, and out-of-plane triggers and for all triggers combined. This is compared to the 
true background, the background from the modified ZYAl method, and the background from the RPF method. (See text for 
details.) The data for all angles relative to the reaction plane have been scaled by 1/3 in the top panel in order to fit on the 
same scale. The fit for the RPF method is to Equationj^to order n = 4 from |A(/>| < 1 and has x^/NDF = 128/90. Bottom: 
Ratios of the background from the RPF and ZYAl methods to the true background. 
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FIG. 13: Top: The true signal for di-hadron correlations in 30-40% Pb + Pb collisions at .^snn = 2.76 TeV for in-plane, mid¬ 
plane, and out-of-plane triggers and for all triggers combined. This is compared to the signal extracted using the background 
from the ZYAl method, the modified ZYAl method, and the background extracted from the RPF method for |A</| < 1 using 
the fit shown in Figure [l^ (See text for details.) The data for all angles relative to the reaction plane have been scaled by 1/3 
in the top panel in order to fit on the same scale. Bottom: Differences between the true signal and the signal extracted using 
the background from the ZYAl method, modified ZYAl method, and the background from the RPF method. 


momenta and we were unable to extend the analysis to 
< 1 GeV/c because in this region the near-side peak 
is too broad even in PYTHIA and distorts the fit. We 
foresee future research using a two-dimensional fit to the 
near-side signal with a Gaussian in A(j) and A 77 in order 
to extend the analysis to lower momenta. 
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FIG. 14: Top: Signal+background for di-hadron correlations in 0-10% Ph -\- Pb collisions at y^SNN = 2.76 TeV in the region 
1.0 < \Ari I < 1.4 for in-plane, mid-plane, and out-of-plane triggers and for all triggers combined. This is compared to the 
true background, the background from the modihed ZYAl method, and the background from the RPF method. (See text for 
details.) The data for all angles relative to the reaction plane have been scaled by 1/3 in the top panel in order to ht on the 
same scale. The ht for the RPF method is to Equationj^to order n = 3 from |A(/>| < 7r/2 and has x^/NDF = 151/140. Bottom: 
Ratios of the background from the RPF and ZYAl methods to the true background. 



FIG. 15: Top: The true signal for di-hadron correlations in 0-10% Pb + Pb collisions at y'sNN = 2.76 TeV for in-plane, mid¬ 
plane, and out-of-plane triggers and for all triggers combined. This is compared to the signal extracted using the background 
from the ZYAl method, the modihed ZYAl method, and the background extracted from the RPF method for \A(j)\ < 'k/2 
using the ht shown in Figure 14 (See text for details.) The data for all angles relative to the reaction plane have been scaled 


by 1/3 in the top panel in order to ht on the same scale. Bottom: Differences between the true signal and the signal extracted 
using the background from the ZYAl method, modihed ZYAl method, and the background from the RPF method. 
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FIG. 16: Top: Signal+background for jet-hadron correlations in 30-40% Ph -F Pb collisions at .^snn = 2.76 TeV in the 

region 1.0 < |A ?7 | < 1.4 for in-plane, mid-plane, and out-of-plane triggers and for all triggers combined. This is compared to 
the true background, the background from the modified ZYAl method, and the background from the RPF method. (See text 
for details.) The data for all angles relative to the reaction plane have been scaled by 1/3 in the tp panel in order to fit on 
the same scale. The fit for the RPF method is to Equation|^to order n = 4 from |A(/| < 7r/2 and has x^/NDF = 193/186. 
Bottom: Ratios of the background from the RPF and ZYAl methods to the true background. 
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FIG. 17: Top: The true signal for jet-hadron correlations in 30-40% Pb + Pb collisions at ^snn = 2.76 TeV for in-plane, mid¬ 
plane, and out-of-plane triggers and for all triggers combined. This is compared to the signal extracted using the background 
from the ZYAl method, the modihed ZYAl method, and the background extracted from the RPF method for \A(j)\ < 'k/2 

(See text for details.) The data for all angles relative to the reaction plane have been scaled 
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using the fit shown in Figure 

by 1/3 in the top panel in order to fit on the same scale. Bottom: Differences between the true signal and the signal extracted 
using the background from the ZYAl method, modified ZYAl method, and the background from the RPF method. 
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TABLE III: Yields Y as defined in Equation scaled by 10“^ from Figure 11 Figure 13 Figure 15 and Figure 17 For the 
true yield the statistical uncertainty is listed first followed by the uncertainty due to the subtraction of the background from 
the underlying event. For the ZYAM method the statistical uncertainty is listed first, followed by the uncertainty due to the 
background level and due to {if and {if. ZYAl uncertainties are propagated assuming 100% correlation between {if and hf and 
no correlation between the uncertainties on the level of the background and on the v„. Uncertainties on the ZYAl method due 
to higher order are not considered but are approximately 10% of the uncertainties due to Df and {if. 


Sample 


near-side Y x 10 

in-plane mid-plane out-of-plane All 

away-side Y x 10 

in-plane mid-plane out-of-plane All 

30-40% 

h-h 

True 

Mod. ZYAl 
Std. ZYAl 
RPF {|A<^| <7r/2) 
RPF (|A<^| <1) 

5.78 ± 0.03 ± 0.13 5.77 ± 0.03 ± 0.14 5.65 ± 0.03 ± 0.13 17.1 ± 0.1 ± 0.2 

6.3 ±5.9 ±1.7 5.7 ±6.0 ±0.3 6.8 ± 6.1 ± 0.9 18.9 ± 4.2 ± 1.2 

4.5 ±2.3 ±1.7 5.5 ±2.3 ±0.3 5.6 ± 2.3 ± 0.9 15.7 ± 1.6 ± 1.2 

5.5 ±0.4 5.7 ±0.3 5.9 ± 0.3 17.0 ± 0.7 

5.7 ± 0.4 5.8 ± 0.4 5.9 ± 0.3 17.4 ± 0.7 

6.74 ± 0.03 ± 0.13 6.72 ± 0.03 ± 0.14 6.52 ± 0.03 ± 0.13 19.9 ± 0.1 ± 0.2 
7.3 ± 5.9 ± 1.7 6.8 ± 6.0 ± 0.3 7.7 ± 6.1 ± 0.9 21.9 ± 4.2 ± 1.2 

5.5 ±2.3 ±1.7 6.5 ± 2.3 ± 0.3 6.5 ± 2.3 ± 0.9 18.7 ± 1.6 ± 1.2 

6.6 ± 0.4 6.8 ±0.3 6.8 ± 0.3 20.1 ± 0.7 

6.8 ± 0.4 6.8 ± 0.4 6.8 ± 0.3 20.4 ± 0.7 

0-10% 

h-h 

True 

Mod. ZYAl 
Std. ZYAl 
RPF (jA<^| <7r/2) 

38.4 ±0.1 ± 0.5 38.2 ±0.1 ± 0.5 37.8 ± 0.1 ± 0.5 114.4 ± 0.2 ± 0.8 

23 ± 26 ± 3 23.2 ± 25.8 ± 0.3 29 ± 26 ± 2 75.5 ± 18.3 ± 0.9 

29.2 ± 9.7 ± 2.8 27.6 ± 9.8 ± 0.3 29.8 ± 9.8 ± 2.3 86.7 ± 7.0 ± 0.9 

35.5 ± 1.2 37.7 ± 0.9 36.0 ± 1.1 109.2 ± 2.5 

44.8 ±0.1 ± 0.5 44.4 ±0.1 ± 0.5 43.8 ± 0.1 ± 0.5 132.8 ± 0.2 ± 0.8 

30 ±26 ± 3 29.9 ± 25.8 ± 0.3 36 ± 26 ± 2 95.7 ± 18.3 ± 0.9 

35.4 ± 9.7 ± 2.8 34.4 ± 9.8 ± 0.3 37.1 ± 9.8 ± 2.3 106.9 ± 7.0 ± 0.9 

41.7 ± 1.2 44.5 ± 0.9 43.3 ± 1.1 129.4 ± 2.6 

30-40% 

jet-h 

True 

Mod. ZYAl 
Std. ZYAl 
RPF (|A<^| <7r/2) 

4.41 ± 0.02 ± 0.10 4.43 ± 0.02 ± 0.10 4.43 ± 0.02 ± 0.10 13.19 ± 0.04 ± 0.17 
4.7 ±5.9 ±1.7 5.2 ±5.9 ±0.3 6.1 ± 6.0 ± 0.9 16.2 ± 4.2 ± 1.2 

3.6 ± 2.2 ± 1.7 4.8 ± 2.3 ± 0.3 4.6 ± 2.3 ± 0.9 13.2 ± 1.6 ± 1.2 

4.3 ± 0.4 4.6 ± 0.3 4.8 ± 0.3 13.5 ± 0.6 

1.67 ± 0.02 ± 0.10 1.67 ± 0.02 ± 0.10 1.70 ± 0.02 ± 0.10 4.96 ± 0.04 ± 0.17 
2.0 ±5.9 ±1.7 2.5 ± 5.9 ± 0.3 3.4 ± 6.0 ± 0.9 8.2 ± 4.2 ± 1.2 

0.9 ± 2.2 ± 1.7 2.2 ± 2.3 ± 0.3 1.9 ± 2.3 ± 0.9 5.2 ± 1.6 ± 1.2 

1.6 ± 0.4 1.9 ±0.3 2.1 ±0.3 5.5 ± 0.6 
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